Pion and kaon elastic form factors in a refined light-front model 
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Within the framework of light-front field theory, we reassess the electromagnetic form factors of 
the pion and kaon. Comparison with experiment is made for the full range of momentum transfer, 
q 2 < 0, including recent data. The light-front model's single regulator mass, tur, of the qq bound- 
state vertex function is initially adjusted to reproduce the weak decay constants, / T and fa, and 
both meson's charge radii, (rv) and (rsr). We study the behavior of these observables under variation 
of the quark masses and find an optimized parameter set, m u — rrid, m s and iur, for which they are 
in sensibly better agreement with experiment than in a previous analysis; a feature also observed 
for the elastic form factors, in particular at small q 2 . This model refinement is important in view of 
an extension to vector and heavy-light mesons. 
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Introduction The light pseudoscalar mesons play a cru- 
cial role in the understanding of low-energy QCD, being 
the lightest strongly bound antiquark-quark states as well 
as the Goldstone bosons associated with chiral symmetry 
breaking. Their static properties have been extensively 
studied — see, e.g., Refs. [THE] — while dynamical features 
have also been investigated theoretically and ex- 

perimentally [30H35j . With respect to the description 
of bound states on the light cone, a detailed review of 
hadronic wave functions in QCD models can be found in 
Ref. [H]. 

In the case of the pion, a great deal has been 
learned experimentally from its electromagnetic form fac- 
tor F 7r (q 2 ) [3"0H3"5] . whereas this is not the case for the 
kaon [371 I3H]- Additional important knowledge about 
the mesons' internal structure can be inferred from their 
valence-quark parton distribution functions 15J. 

The theoretical framework we adopt is the light-front 
field theory formalism [3BJ. More specifically, we here 
ameliorate the light-front approach first introduced in 
Refs. [U [TJJ , where two classes of qq bound-state mod- 
els for the Bethe-Salpeter amplitude of the pion must be 
distinguished: the nonsymmetric |17j and the symmet- 
ric |18j vertex model. The light-front component J + of 
the electromagnetic current has been successfully used to 
calculate elastic form factors [TJ [TBI I39H43] . For the non- 
symmetric irqq and Kqq vertex models, the components 
of the current are conveniently obtained in the Drell-Yan 
frame. We recall that on the light cone the bound state 
wavefunctions are defined on the hypersurface x°+x 3 = 
and are covariant under kinematical boosts due to the 
stability of Fock-state decomposition [TT] . 

In this Brief Report, we solely consider the nonsym- 
metrical qq bound-state vertex function with the inten- 
tion to optimize and unify the parameter set which simul- 
taneously reproduces the pion and kaon decay constants, 
charge radii and their electromagnetic form factors. For 
the latter, our numerical results are compared with ex- 
perimental data up to 10 GeV 2 in order to explore the 
validity of the model at large q 2 transfer. 



The model We briefly summarize the light-front model 
including the Ansatz for the nonsymmetric vertex func- 
tion for the pseudoscalar bound states. The covariant 
electromagnetic form factor is defined as (suppressing 
color and flavor indices), 

(p + P '),F Mo _{q 2 ) = (Mo-ip'^Ma-ip)) , (1) 

in which Mq- = ir, K denotes the light pseudoscalar 
mesons, q = p—p' and J M = e q ipq'y^ipq is the electromag- 
netic current. In the impulse approximation, the form 
factor is given by a triangle diagram which represents an 
amplitude via a single integral: 



(Pk+p'JFm (q 2 ) =e q N c 



T(k,p')T(k,p) 



■H'-«o-« ' "1'--J ( 27r )4 

x Tv\s q (k) l5 S q (k - p')^S q (k - p) 75 +[q^q], (2) 



where S q 1 (p) = j> - 



te, q = u, q = d, s, denotes the 



quark propagator with quark masses m q , e q is the quark's 
charge and N c = 3 is the number of colors. The bracket 
[q <H- q] is a shorthand for the exchange of the quark and 
antiquark in the integral. In the following, we work in 
the Breit frame with the Drell-Yan condition, such that 
= bo, 9/2, 0,0), p'^ = (p' ,-q/2,0,0), p Q = p' and 
?p = (0,g,0,0). 

After transformation to light-cone variables, using the 
plus component, J + , of the electromagnetic current, the 
first integration is over kT (the null-plane energy) whose 
pole contribution is fc _ = (k± + m q — ie)/k + . It was 
shown that k + is limited to two integration intervals: the 
valence contribution, constrained by < k + < p + , while 
the nonvalence contribution is restricted to p + < k + < 
p' + [T71 118j . Thus, on the light front, two distinct terms 
contribute to the electromagnetic form factor: 



(3) 



where F]^ 



(q 2 ) and Fj^ _ (q 2 ) are the valence (qq) and 
nonvalence (pair production) contributions, respectively. 
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Furthermore, it can be shown that for a non symmet- 
ric bound-state vertex function, T(k,p), asymmetric un- 
der momentum exchange of the quark and antiquark, the 
second interval does not contribute to the electromag- 
netic form factor [17] . Within a quark-meson interaction 
model with effective coupling, 



£ 



7T • q~/ 5 Tq 



a nonsymmetric vertex function is 



T(k,p) = 



N 



(p-ky 



(4) 



(5) 



where, m — (m q + rriq)/2, /m _ is the weak decay con- 
stant, rriR is the regulator mass and N is an overall 
normalization. Thus, after kr integration and change 
of variables, x = k + /p + , where x is the momentum 
fraction carried by the quark in the infinite-momentum 
frame, the elastic form factor is given by the integral 



F + - 

r M, — 



e q N c 



d 2 k±dx 



fej.)*(x, k±)9(l - x)9(x) 



2p+ J (2vr) 3 x 
Tr[fJ + ™ g -b 5 (^ ~i>' + m^+W -i> + m g ) 7 5 
+ [q++q] . 



k-=k- 

(6) 



The subscript "+" on the form factor in Eq. (6| is a 
reminder that we employ the J + component of the elec- 



tromagnetic current. For M - = 
are q = u and q = d, whereas for M 



the quark flavors 
K + we have 



q = u and q — s and [q ■<-> q] is again an abbreviation 
for both quark and antiquark contributions to the elastic 
form factor. Moreover, after k~ integration, a light-front 
wave function emerges which for a nonsymmetric qq ver- 
tex function is defined as: 



*(x, kj 



N 



f Mo _ {l-xy{m? M -Ml){r 



Ml) 



where Mr — M(m 2 ,m R ) is given by, 
\ , (p-k)l + 



Mi 



kl 



(l-x) 



Pi 



(7) 



(8) 



Mq = M 2 (m 2 ,m 2 ) is the free mass operator, tum _ is 
the pseudoscalar meson mass and the normalization con- 
stant N obeys the condition F^ f _ (0) = 1. 

In addition, we calculate the pseudoscalar weak decay 
constant with the same bound-state vertex model intro- 
duced in Eq. 



(0\A tl (p)\M -(p))=%f Mo _p M 



(9) 



where — i^q^^^q- After Dirac algebra and k 
integration, the pseudoscalar decay constant reads, 



N, 



d 2 k±dx 
(2n) 3 x 



[4xmq+4m q (l-x)]^(x, k±). (10) 



We remind that the charge radius is determined via, 



dF Mo _ (g 2 ) 
dq 2 



(11) 



q 2 =0 



Numerical Results We have three model parameters: 
the regulator mass, rriR, and the quark masses, m u = m c i 
and m s . The main aim of this work is to jointly analyze 
the pion's and kaon's elastic form factors, decay constants 
and charge radii in order to determine more accurately 
the model's quark masses in view of future applications 
and to test whether a single mass scale, mn, can satisfac- 
torily describe experimental data for both light mesons. 

We initially consider the parameters of Ref. [55] > i.e. 
m u = m d = 220 MeV, m s = 419 MeV, = 140 MeV, 
m K+ = 494 MeV and m R = 0.946 GeV, in which the 
elastic form factor ratio of the kaon and pion were cal- 
culated. This parameter set describes rather well mea- 
sured observables for the pion but less so for the kaon: 
(r 7r +) = v / ( 7 ' 2 + ) = 0.672 fm, which coincides with the 
experimentarvalue, (r^+ p ) = 0.672 ± 0.008 fm [M], and 
U = 101 MeV (/r p 7V 2 = 92 -21 ± 0.14 MeV @7]); 
whereas (r K +) = yJ{r 2 K+ ) = 0.71 fm ((r^) = 0.560 ± 
0.031 fm [37]) and f K = 129 MeV (f^'/V 2 = H0.4± 
0.6 MeV 07J) indicate that the model is not well adjusted 
for mesons with strangeness content. 

We vary the masses of the constituent quarks and si- 
multaneously tune mn toward a common value for both 
the pion and the kaon. Clearly, the regulator mass acts 
as a cut-off in the triangle diagram but also defines a 
physical length scale. Therefore, one would not expect 
m R to be equal for both mesons. However, we here in- 
sist on having a minimal number of parameters which 
still yields a satisfying reproduction of all available data. 
Once the light quark masses in the model are fixed, we 
can consider heavier mesons as well as 1~ mesons [48]l49] 
and compute their electromagnetic properties [50] , 

In Fig. [I] we plot the charge radii (r„) and {rx) as 
a function of m u — m^, where m R — 1.0 GeV and 
m s = 510 MeV are kept fixed, whereas in Fig. [2] the 
functional behavior of the charge radii in dependence of 
rriR is shown with m u = 220 MeV and m s = 510 MeV 
fixed. The dependence is in both cases nonlinear and 
somewhat more pronounced for variations of the quark 
mass than of the regulator mass. We remark that a larger 
regulator mass results in a smaller charge radius, as ex- 
pected. The strange constituent quark mass has been 
readjusted from its value m 8 = 451 MeV [33 to obtain 
a better agreement with the kaon charge radius while 
keeping mn — 1.0 GeV. The value m s = 510 MeV we 
here choose is in agreement with a definition of the Eu- 
clidean constituent quark mass derived from solutions of 
Dyson-Schwinger equations |51[ 152]. 

Similarly, in Fig. [3] the decay constants f n and fx are 
shown as a function of the quark mass (to_r = 1.0 GeV) 
and in Fig. [4] as a function of the regulator mass (where 
m Utd = 220 MeV and m s = 510 MeV). We note that 
the decay constants are also more sensitive to varia- 
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FIG. 1: Pion and kaon charge radii, v/ (r 2 ), as a function of 
m u = rrid\ m a = 0.51 GeV and mj = 1.0 GeV are fixed. 
Upper and lower dotted lines mark experimental values for 
the pion and kaon, respectively. 




m R [GeV] 

FIG. 2: Charge radii, (r 2 ), of the pion and kaon as function 
of tur with m u — md = 0.22 GeV and m s — 0.51 GeV fixed. 
Horizontal dotted lines as in Fig. [l] 

tions of m q than of m,R, excepting in the case of the 
strange quark where a departure of about 10% from 
rn s = 510 MeV does not significantly alter Jk ■ For the 
pion, m R — 1.0 GeV yields the best adjustment to the 
experimental values of the decay constant and charge ra- 
dius; since it also occurs to best describe the pion's elas- 
tic form factor in Figs. [5] and [7J as discussed below, we 
definitely set tjir = 1.0 GeV and use this value in calcu- 
lations of the kaon's properties. Moreover, while (rv) and 
(rji) decrease with increasing quark mass, the opposite is 
true for /„. and fx , which is an expression of the Tarrach 
relation, (r M[) _ ) ~ l/m q and f Mg _ ~ l/{r Mo _) [S3]. 

We thus require m R = 1.0 GeV, m u = m d = 220 MeV 
and m s = 0.51 GeV as our reference values and ob- 
serve their implications for both pseudoscalars' elastic 
form factors. In case of the pion, this mass choice in the 
light-front model reproduces very well the experimental 
data [SHI [3TJ |3"3"H3"5] , while the lack of data on the kaon's 
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FIG. 3: Weak decay constants: f n as function of m u = mj 
(solid line) and fx as function of m u (dashed line with m s = 
0.51 GeV) and m a (dot-dashed line with m u — 0.22 GeV); 
m,R = 1.0 GeV in all cases. The upper and lower dotted lines 
denote experimental values for fx and /„■, respectively. 
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FIG. 4: The weak decay constants, f n and fx, as a function 
of the regulator mass, ttir, with m u = md = 0.22 GeV and 
m s = 0.51 GeV fixed. Horizontal dotted lines as in Fig. [3j 

elastic form factor and associated large error bars do not 
allow for a satisfactory comparison. 

For the sake of completeness, we illustrate the strong 
sensitivity of the elastic form factors to tur. As can be 
seen from Figs. [5] and [6j acceptable values with respect 
to the experimental data on F^(q 2 ) and F^(q 2 ) lie in the 
interval, 0.8 GeV < mji < 1.3 GeV, which coincides with 
our privileged value mu = 1.0 GeV. 

Next, the model dependence on the constituent masses 
are explored in Figs. [7] and [8j For the pion, we observe 
that this dependence is asymmetric and more strongly 
pronounced for smaller quark-mass values. For instance, 
compared with the reference mass m u — 220 MeV, a 
70 MeV lighter constituent quark yields a soft elastic 
form factor that strikingly deviates from the experi- 
mental data, whereas the same calculation with m u = 
300 MeV results in just a slightly harder form factor, in 
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FIG. 5: Elastic form factor of the pion as function of ran and 
with m u = mj = 0.22 GeV. Experimental data: Refs. [301131] 
(solid triangles), Ref. [33] (solid diamonds), Ref. [34] (solid 
circles) and Ref. [35] (solid squares). 



FIG. 7: Elastic form factor of the pion as function of m Ut d = 
0.22 GeV and mu = 1.0 GeV. Experimental data as in Fig. [5] 
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FIG. 6: Elastic form factor of the kaon as function of ran 
with m u = rud = 0.22 GeV and m^ = 0.51 GeV; experimental 
data from Refs. [371155] . 

particular for Q 2 = —q 2 < 0.2 GeV 2 . This is not sur- 
prising, as this momentum range probes the static fea- 
tures of the mesons where a constituent quark mass of 
200 — 350 McV is appropriate [53] • For larger values of 
q 2 , however, this is not the case and a simple model will 
fail. This asymmetric behavior with respect to the light- 
quark mass is not observed in Fig. [HJ where the kaon's 
dynamical features are dominated by the strange mass. 

These observations on the elastic form factors confirm 
our choice of mass parameters, m,R = 1.0 GeV, m u = 
nid = 220 MeV and m s = 510 MeV, for which we obtain: 
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Conclusions We reassessed the light-front model [TT] 
ITS] in view of recent data on F 7T (q 2 ) and due to the need 
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FIG. 8: Elastic form factor of the kaon in dependence on the 
quark mass, m u , where ran = 1.0 GeV and ra s - = 0.508 GeV 
are fixed; experimental data are from Refs. [371 138] . 



of a better determination and restriction of the model 
parameters. The regulator mass, ttir = 1-0 GeV, us- 
ing a nonsymmetric vertex model for the bound state, 
is found to simultaneously satisfy the experimental data 
on the space-like elastic form factors, the weak decay 
constants and the charge radii of the pion and the kaon 
within reasonable theoretical uncertainties. The numeri- 
cal results show that the model significantly breaks down 
for < 0.8 GeV and for m,R > 1.3 GeV. We have also 
studied the model dependence on the constituent masses, 
m u = rrid and m s , showing that Tarrach's relation is sat- 
isfied and confirming the range of mass values commonly 
chosen within the light-front model. These parameter 
values are useful for, e.g. heavy-meson and vector decay 
constants [SD] or heavy-to-light transition form factors. 
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